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We discuss the properties of nonlinear localized modes in sawtooth lattices, in the framework of 
a discrete nonlinear Schrodinger model with general on-site nonlinearity. Analytic conditions for 
existence of exact compact three-site solutions are obtained, and explicitly illustrated for the cases 
of power-law (cubic) and saturable nonlinearities. These nonlinear compact modes appear as contin¬ 
uations of linear compact modes belonging to a flat dispersion band. While for the linear system a 
compact mode exists only for one specific ratio of the two different coupling constants, nonlinearity 
may lead to compactification of otherwise non-compact localized modes for a range of coupling ra¬ 
tios, at some specific power. For saturable lattices, the compactification power can be tuned by also 
varying the nonlinear parameter. Introducing different on-site energies and anisotropic couplings 
yield further possibilities for compactness tuning. The properties of strongly localized modes are 
investigated numerically for cubic and saturable nonlinearities, and in particular their stability over 
large parameter regimes is shown. Since the linear flat band is isolated, its compact modes may be 
continued into compact nonlinear modes both for focusing and defocusing nonlinearities. Results 
are discussed in relation to recent realizations of sawtooth photonic lattices. 

PACS numbers: 63.20.Pw, 42.82.Et, 78.67.Pt, 03.75.Lm 


I. INTRODUCTION 

The creation and manipulation of nonlinear local¬ 
ized lattice excitations (“discrete solitons” or “discrete 
breathers”) is of interest in many areas of physics [1], 
and in particular within nonlinear optics [2]. In many 
contexts, one strives for achieving optimal spatial local¬ 
ization. In the presence of linear dispersion (diffraction), 
resulting, e.g., from a linear (harmonic) coupling between 
neighboring sites, the localization will generically be ex¬ 
ponential. However, non-generic situations may appear 
in presence of nonlinear dispersion, e.g., when the inter¬ 
site interaction also has a nonlinear contribution. In such 
situations, a proper control of the nonlinear effects may 
lead to a complete cancellation of the linear coupling be¬ 
tween sites, resulting in exact “discrete compactons” fully 
localized at a small number of sites, with strictly zero tail 
(for examples, see, e.g., [3-5]). The experimental realiza¬ 
tion of these specific nonlinear interactions poses how¬ 
ever serious challenges, and so far no direct observation 
of compactification of nonlinear lattice modes has, to the 
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best of our knowledge, been reported. 

On the other hand, it is well known that in certain lat¬ 
tices geometries, frustration effects make possible the ex¬ 
istence of strictly compact modes even in absence of any 
nonlinearity. Such modes can be seen to result from de¬ 
structive interference of couplings between different sites, 
and are typically associated with a flat linear dispersion 
band. A number of different examples of this scenario 
are discussed, e.g., in Refs. [6-9], and references therein. 
In particular, recent experiments have being performed 
in two-dimensional (2D) Lieb photonic lattices, where 
different transport regimes [10] as well as the excitation 
of the localized flat band mode [II, 12] were observed. 
Maybe the most well known example is the 2D Kagome 
lattice, which supports modes strictly localized on six 
sites along a hexagonal geometry (see, e.g., [13]). As was 
shown in Ref. [14], these compact modes exist also in 
presence of a pure on-site nonlinearity, and indeed con¬ 
stitute an effective ground state of the system in the 
small-power regime for a defocusing Kerr (cubic) non¬ 
linearity. However, above a rather small threshold power 
the compact modes destabilize, and transform into stable 
standard localized modes with exponential tails. 

The compactification of linear modes due to competing 
linear interactions may appear also in one-dimensional 
(ID) settings, and a number of examples have been illus¬ 
trated, e.g., in Refs. [6, 8]. One of the most studied ID 
lattices is the sawtooth chain; see, e.g., [6, 7, 15-17] (also 
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called A chain [18] or triangle lattice [ 8 , 19]). Due to 
its construction from corner-sharing triangles, it is often 
considered as a ID analogue of the Kagome lattice, and is 
one of the simplest possible lattices allowing for compact 
linear modes and a flat band, which makes it particu¬ 
larly attracting for experimental realizations. A compact 
mode in a similar setting was used to experimentally ob¬ 
serve a Fano resonance [20]. Very recently, the linear 
properties of a sawtooth photonic lattice, created using 
the femtosecond laser written technique, were also stud¬ 
ied experimentally [21]. Furthermore, the dependency of 
coupling constants of these lattices on polarization [ 22 ] 
provides ample and tunable opportunities beyond geo¬ 
metric considerations. This is important since, in con¬ 
trast to the Kagome lattice, compactification in a linear 
sawtooth lattice appears only at one specific (irrational) 
ratio between the two different coupling constants (hori¬ 
zontal and diagonal), which again poses an experimental 
challenge for fine-tuning. However, as we will show in 
this work, nonlinearity allows the existence of compact 
modes for different coupling ratios, which certainly would 
facilitate the observation of highly localized states in a 
real experiment. 


The aim of the present paper is to illustrate how non- 
linearities and other effects can be used for the purpose 
of facilitating the compactification tuning of localized 
modes in sawtooth lattices. In Sec. II we describe the 
general model, using the discrete nonlinear Schrbdinger 
(DNLS) equation with a general on-site nonlinearity and 
linear coupling constants representing the sawtooth ge¬ 
ometry. We discuss the linear dispersion relation and its 
flat-band structure, and obtain the general conditions on 
the geometry for existence of 3-site compactons as exact 
stationary solutions. In Sec. Ill we specialize to power- 
law nonlinearities, and obtain an expression for the power 
needed for compactification, valid for a whole range of ra¬ 
tios for the linear coupling constants. For the physically 
most interesting case of cubic Kerr-nonlinearity (focusing 
or defocusing), we perform a detailed numerical anal¬ 
ysis of localized modes (not restricted to compactons) 
in Sec. Ill A, illustrating explicitly how certain, linearly 
stable, localized modes compactify when the appropri¬ 
ate relation between coupling ratios and power is ful¬ 
filled. Another physically important case of a saturable 
nonlinearity is analyzed in Sec. IV. The saturability in¬ 
troduces another parameter which can be employed, to¬ 
gether with the power, for compactification tuning, and 
we illustrate the scenarios in different regimes analyt¬ 
ically as well as numerically, particularly showing the 
existence of strongly localized, stable modes over large 
parameter regimes. Finally, Sec. V discusses two addi¬ 
tional generalizations of the sawtooth lattice with cubic 
nonlinearity, yielding further possibilities for compactifi¬ 
cation tuning: on-site energy alternations (Sec. V A) and 
coupling anisotropies (Sec. VB). Conclusions are given 
in Sec. VI. 



FIG. 1: (Color online) Geometry of the sawtooth lattice with 
its compact mode (white circles imply zero amplitudes). 


II. MODEL AND GENERAL COMPACTON 
CONDITIONS 

We consider the following form of the discrete nonlin¬ 
ear Schrbdinger equation on a sawtooth lattice sketched 
in Fig. 1, 


iUn T ^ “t” T/(jt^n j )tin — 0 i (1) 

m^n 

where n = 1,2,...,V. The coupling Jn,m denotes the 
interactions between the site n and its connected lattice 
sites m, with diagonal (“tip”) couplings J* and horizon¬ 
tal (“line”) couplings J;. The function f{x), describing 
a general on-site nonlinearity with strength 7 , is at this 
point arbitrary. We may put /(O) = 0, since any ho¬ 
mogeneous linear on-site potential may be removed by 
a gauge transformation (extensions to non-homogeneous 
on-site potentials will be discussed in Sec. V). Generally, 
the parameters J/, Jt and 7 may be complex-valued (de¬ 
scribing effects of gain and/or loss), but in our explicit 
examples in Secs. III-V we will only consider real-valued 
cases. We should remark that the geometry in Fig. 1 
apparently differs from that usually associated with saw¬ 
tooth lattices [15, 18]. Instead of taking equally oriented 
triangles, we take the triangle tips to point alternately up 
and down from the horizontally line-coupled sites. In this 
way, we avoid the otherwise physically unavoidable inter¬ 
action between neighboring vertices, describing a more 
appropriate experimental realization of a photonic saw¬ 
tooth lattice [21]. This obviously does not change the 
mathematical model, which in both cases only contains 
the two nearest-neighbour couplings Jt and J/. 

Looking for stationary solutions of the form 

Unit) = anC"^* , (2) 

the first term in Eq. (1) is replaced with —Xun- Assuming 
the line sites to be of odd n, whereas for the tips n is even, 
then leads to the following set of equations in the linear 
regime (qlo^l^ = 0 ): 


. j Jiian +2 + an- 2 ) forn odd , 
Aa„ — Jt(a„-|-i+a„_i)-|-< 

(J tor n even . 

The linear dispersion relation [ 6 , 7, 21] is computed by 
assuming a propagating plane wave a„ = Ae*^", obtain¬ 
ing for the linear frequencies, rescaled with the horizontal 
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FIG. 2: (Color online) Dispersion relation (3) for the sawtooth 
lattice for different coupling ratios: J= 0.1 (0.5, 1, \/2, 2, 2.5) 
in thin black (thick black, thin orange, thick red, thin dashed, 
thick dashed). 

coupling constant J;, 

A±(fc;J) = A±M= (3) 

cos(2fc) ± ^Jcos^{2k) + 2,P [1 + cos(2fc)] . 

Here we have also introduced the rescaled, dimensionless 
coupling ratio J = Jt/Ji, which is the essential system 
parameter that determines its qualitative physical prop¬ 
erties in the linear regime. 

As illustrated in Fig. 2, the lower band becomes flat at 
the particular coupling ratio J = ^/2 for which A_ = —2. 
It is important to note that the flat band existing at 
these parameter values is isolated from the upper band 
(A+ = 4cos^ k) with a bandgap of 2. Thus, in contrast 
to the Kagome lattice studied in Ref. [14], where the 
flat band is always attached to the edge of a dispersive 
band, we may here expect a resonance-free continuation 
of compact modes belonging to the flat linear band for 
focusing as well as for defocusing nonlinearities. 

To obtain a condition for existence of compact station¬ 
ary solutions (2) of Eq. (1) for general nonlinearities /(cc), 
we write an ansatz for the solution as a 3-site compacton 
(analogously to the linear flat-band modes [6, 7, 21]), 
with a„p = A for the central, 4-bonded, “line site” (red 
disk in Fig. 1), and a„(,±i = aA for the two neighboring, 
2-bonded, “tip sites” (blue in Fig. 1). The condition to 
have zero amplitude on the neighboring line sites ng ± 2, 
and thus also for the rest of the lattice, then gives imme¬ 
diately an expression for the amplitude ratio a in terms 
of the ratio between the coupling constants: 

a = -1/J . (4) 

For the central line site ng, Eqs. (l)-(2) give the following 
expression for the rescaled frequency A/J;: 

AM = -2+^/(lAl2) . (5) 

•Ji 


The corresponding equations for the two tip sites ng ± 1 
yield 

AM = -J2 + ^/(1H/J12) . (6) 

•Jl 

Eliminating A from Eqs. (5)-(6) gives an implicit, general 
condition to determine the central-site intensity in 

terms of the coupling ratio J and the ratio y/J; between 
nonlinearity strength and horizontal coupling constant: 

T^f(\A\^)- f{\A/J\^))=2-J^ . (7) 

Obviously, in the linear limit yJA]^ —>■ 0 and the left-hand 
side of Eq. (7) vanishes. This implies that Jt —>■ ±y/2Ji, 
recovering the linear flat-band condition. Using Eq. (4), 
we may also determine the total power, P = 
of the compacton as 

P = {1 + 2\1/J\'^)\A\^ . (8) 

Apart form the power, the other conserved quantity of 
DNLS-like systems is the Hamiltonian 

H=Y1 ( 9 ) 

n^m n 

for the canonical variables The function 

F{\un\‘^) defines the nonlinearity of the system {f{x) = 
F'{x)) and will be specified later on. (Note that we 
chose here a sign-convention in accordance with the sign- 
convention used in the definition of the frequency A in (2). 

As illustrated below, stable compactons may then appear 
as ground states for defocusing nonlinearities.) 

In order to compare how compactons are related to 
eventual ordinary localized stationary modes, which have 
nonvanishing tails, we implement a Newton-Raphson it¬ 
erative method to find numerical solutions at a given ac¬ 
curacy. We will look for fundamental solutions centered 
at tip or line sites or in between sites, using periodic 
boundary conditions to avoid surface effects. A typical 
lattice size used numerically is A^ = 40, while also check¬ 
ing the consistency of the results by varying the chain 
length. We generally use J/ = 1 in all numerical calcula¬ 
tions. 

For any nonlinear stationary mode {an} fulfilling (1)- 
(2), we perform a standard linear stability analysis [23]: 
we use the ansatz Un{t) = [a„ -I- e„(t)]e*''*, and linearize 
equations in e„(t), giving access to the spectrum {wi} 
of perturbations e*“‘‘. We define a stability parameter 
g = Maxjlm(wi)j, which will be nonzero for unstable 
modes, whereas 5 = 0 for stable solutions. To compare 
the localization properties of different modes, we use the 
participation number defined as 

1 for a single excited site, 

N for extended homogeneous profiles. 

This parameter gives a measure of the effective size of a 
given profile. In particular, for the compacton solution 


R = 


p2 


E 


N I 
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this formula reduces to 


Rc 


\J\^ + 2 

W+2 


2 for J ^ 0 , 

1 for J —>■ oo . 


( 10 ) 


From this, it also follows that the maximum compacton 
participation number i?c = 3 would appear for J —> 1, 
when the profile would be equally distributed between all 
3 sites. For the linear case, Rc{\J\ = V^) = 8/3. 


III. POWER-LAW NONLINEARITIES 


Specializing first to a general, pure power-law nonlin¬ 
earity, we have f{x) = x'^. Plugging this into Eq. (7) 
yields that an exact 3-site compacton exists when the 
following relation, between the central-site intensity and 
the coupling constants, is fulfilled: 


2 -J 2 

Ji ~ 1-1/kP- ■ 


( 11 ) 


Using Eq. (8), we find that the compacton condition can 
alternatively be written in terms of the total power as: 


P=(2+|J|2) 


U(i^P"-i); 


( 12 ) 


For power-law nonlinearities, F{x) = /{a -\- 1) in 

Hamiltonian (9); for the compact mode we find 


Hc=y= (13) 

^JP 2 \ / |j|2p x(-+i)/2 

2 + I J|2 ^ a + 1 'v / V2 + \j? ) 


A. Localized modes for cubic nonlinearities 


For the special case cr = 1, which corresponds to a 
standard cubic (Kerr) nonlinearity, and real-valued cou¬ 
pling constants, Eq. (12) considerably simplifies. We can 
define an effective nonlinearity parameter as 


Py _ 4-J4 

~ Ji - 1 ■ 


(14) 


For Kerr nonlinearities, F is the only free parameter. 
Therefore, for focusing nonlinearities F > 0, whereas 
F < 0 for the defocusing case. We can also calculate the 
rescaled frequency of the (nonlinear) compacton, given 
by 


T Ac J4 _ 2 


(15) 


The behavior of both F and Ac is shown in Fig. 3 with 
a solid black and dashed line, respectively. Furthermore, 



FIG. 3: (Color online) Properties of the compacton: effective 
nonlinearity F (solid), rescaled frequency Ac (dashed) and sta¬ 
bility parameter g (dotted) vs. J. The linear band extension 
is shown with shaded areas. The analytical expressions for F 
and Ac are equivalent to those of a three-core coupler stud¬ 
ied in [24] for different regimes: 0 < J < 1 “branch c”; 
1 < J < e-)- “branch b”; J > %/2 “branch d”. 


we show the stability parameter g (red dotted) vs. J. 
Whenever g > 0, the solution is unstable. The lin¬ 
ear band extensions are plotted as shaded background. 
When the linear compacton condition J — \p2 is fulfilled, 
obviously F(J) = 0 and Ac = —2. At this value, we ob¬ 
serve that the lower band width becomes strictly zero, 
corresponding to the perfectly flat band of this lattice, 
occurring for this special geometric relation. 

For J > -\/2, the linear band broadens but a compact 
mode still exists for negative F, with a frequency below 
the lower band. As is seen, the mode stays linearly stable, 
which was also confirmed by direct numerical integration 
of Eq. (1). In the large-J limit, |a| decreases and the 
compacton, having its dominating amplitude on a line 
site, approaches a highly nonlinear single-site localized 
mode of an ordinary defocusing ID chain (as the line- 
couplings become negligible next-nearest neighbor cou¬ 
plings), which is always stable. 

On the other hand, for decreasing J < \f2 compactons 
instead occur for a focusing nonlinearity (F > 0) and 
with the frequency entering into the gap above the lower 
band. As can be seen, the mode initially remains linearly 
stable but becomes unstable in a bifurcation at J ~ 1.27. 
The origin of this instability can be understood by com¬ 
paring with the well-known properties of the stationary 
solutions of a three-core coupler with focusing nonlin¬ 
earity (with open boundary conditions, i.e., without any 
surrounding lattice), studied, e.g., in Ref. [24]. The an¬ 
alytical forms of the compactons in the sawtooth lattice 
are mathematically equivalent to those of the symmetric 
three-core eigenmodes given in [24]; however, the stabil¬ 
ity properties may differ due to additional possible reso¬ 
nances with the surrounding lattice. As the correspond¬ 
ing modes in the focusing three-core coupler (“branch b” 
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with the notation of [24]) are always stable, the observed 
instabilities for the compacton in the focusing sawtooth 
lattice must originate in resonances between internal os¬ 
cillations of the compacton core and linear oscillations 
of the surrounding chain. Decreasing J further towards 
the singularity at J = 1, stronger instabilities appear as 
the compacton frequency becomes positive and enters the 
upper linear band. Compactons with positive frequency, 
in a strongly nonlinear focusing lattice, are found to be 
always unstable. 

From the above, we reach the important conclusion 
that there are stable compact nonlinear modes in direct 
continuation of the linear flat-band modes for both fo¬ 
cusing and defocusing nonlinearities. At the singularity 
J = 1 both r and Ac change sign, and we see that for 
J < 1, compact modes with negative effective nonlinear¬ 
ity and frequency are always unstable. Noting that the 
anti-phased compact modes in this defocusing regime are 
mathematically equivalent to three-core in-phase modes 
with focusing nonlinearity, these instabilities are equiv¬ 
alent to those described in Ref. [24] for the solution 
“branch c” (corresponding to main amplitudes at the two 
tip sites). Thus, this instability originates from an inter¬ 
nal resonance between oscillation modes of the three-site 
compacton core. So to conclude this part, we see that 
for any sawtooth lattice geometry we can find a compact 
mode centered on a line site. This mode is stable for a 
wide window in parameter space but destabilizes when 
the surrounding tip site amplitudes become comparable 
to, or larger than, the amplitude of the central line site. 

The nonlinear compactons described above are not iso¬ 
lated solutions for a chain with a given J, but rather be¬ 
long to certain families of, generally non-compact, non¬ 
linear localized modes which may be tuned to compact¬ 
ness by varying the effective nonlinearity (e.g., by power 
tuning) so to fulfill the compacton condition (14). We 
first illustrate this scenario for an effectively defocusing 
nonlinearity, T < 0, in Fig. 4. The relations of A, R, 
g and H vs. J, for four families of ordinary localized 
modes with F = —2.5, in comparison to the compact 
mode with its respective F given by (14), are shown in 
Fig. 4 (a)-(d). Here, the thick (thin, dashed, dotdashed) 
line represents the compact (odd line, odd tip and even 
line) mode; some examples for these modes are shown in 
Figs. 4 (e)-(g), with parameter values corresponding to 
the respective symbols. The odd line mode is centered at 
a line site and has the same symmetry as the compacton, 
while the even line mode, as well as the odd tip mode, 
have different symmetries being centered around or at 
a tip site (the even line mode has its main amplitudes 
at the surrounding line sites, while the odd tip mode is 
primarily localized at the central tip site). As usual in 
cubic lattices, the even line mode is unstable everywhere. 
An example for its profile a„ is shown in Fig. 4 (e). We 
chose parameter values, where its participation number 
R is very close to the compact mode, but, there is a 
difference in frequency and Hamiltonian, and the profile 
shows broader exponential tails. The same happens for 
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FIG. 4: (Color online) Properties of defocusing nonlinear 
modes: (a) frequency A, (b) R, (c) stability parameter g and 
(d) H vs. J for the compact, odd line, odd tip, even line and 
intermediate modes in thick, thin, dashed, dotdashed and dot¬ 
ted lines, respectively. The linear band in (a) is shown with 
a shaded area. In the inset of (d) we show the difference be¬ 
tween H for the odd line and tip mode. For all non-compact 
modes F = —2.5. The symbols represent the parameter values 
of the exemplary modes shown below: (e) an of the even line 
mode, (f) fln of the odd tip mode, (g) a„ of the odd line mode, 
(h) On of the intermediate mode. Red and blue filled circles 
represent the amplitude at line and tip sites, respectively. 


the odd tip mode Fig. 4 (f), which is also not directly 
connected to the compact mode. Finally, the family of 
odd line modes possesses one compact mode for the spe¬ 
cific parameter relation (F, Ac, T), which is shown in Fig. 
4 (g). Thus, the nonlinear compact modes appear as 
special cases of the continuous family of stable, exponen¬ 
tially localized, modes centered on a line site, which can 
be tuned to compactification at these specific parameter 
values. 

For larger J, the odd line mode has an interval of in¬ 
stability around 2.43 ^ J ^ 5, where it changes its form 
(upper bifurcation not shown in Fig. 4(c)). The tip sites 
next to the central line site grow and pass its amplitude. 
A similar instability interval appears for the family of 
odd tip modes for smaller J, 1.7 ^ J ^ 2.05. In the in¬ 
terval between the two bifurcation points where the odd 
tip modes stabilize and the odd line modes destabilize, 
a family of intermediate, asymmetric modes appears (cf. 
Fig. 4(h)). This unstable mode, shown in Fig. 4 with a 
dotted line in its regime of existence, is responsible for 
transferring the instability between the odd tip and line 
modes, and merges with the odd modes at the bifurca- 
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FIG. 5: (Color online) Properties of focusing nonlinear modes: 
(a) frequency A, (b) R, (c) stability parameter g and (d) H 
vs. J for the compact mode (thick black), the compactifying 
odd line gap mode when P = 1.2 (thin red), and the non- 
compactifying odd line mode above the spectrum when F = 5 
(dashed green). Examples of stable modes are shown below: 
gap modes in (e)-(f) and modes above the spectrum in (g)-(h). 


tion points. In the inset of Fig. 4(d) we show the dif¬ 
ference between H for the odd line and tip mode, where 
the crossing hints at this stability exchange scenario and 
the existence of the intermediate solution. Note however, 
that this regime is rather far from the compactification 
regime. 

The analogous compactification of a nonlinear local¬ 
ized, stable gap mode for focusing nonlinearity is illus¬ 
trated in Fig. 5 (for simplicity, we include in this figure 
only modes with main peak at a single line site). It is 
important to stress that the compactifying mode, with 
its sign changing amplitudes (Fig. 5 (e), (f)) is not iden¬ 
tical to the fundamental odd line mode with amplitudes 
of the same sign (Fig. 5 (g), (h)) and frequency above the 
spectrum, and they may both exist as simultaneously sta¬ 
ble solutions. The condition to have the compactifying 
mode stable at the point where it compactifies requires 
the nonlinearity to be relatively weak, F < 2.3 (cf. (14)). 
On the other hand, the fundamental odd line mode ex¬ 
ists as a stable solution also for larger F, as illustrated 
in Fig. 5 for F = 5, where the compact mode also has 
a frequency above the linear spectrum but is unstable. 
The properties of the compactifying modes for stronger 
focusing nonlinearities will be discussed further for the 
saturable nonlinearities in Sec. IV (cf. Fig. 8 ). 


IV. SATURABLE NONLINEARITY 


Here, we consider a saturable nonlinearity, taken in 
equivalent form as in, e.g.. Ref. [25] and references 
therein. 


Thus, in the small-amplitude limit, the Kerr nonlinear¬ 
ity of Sec. Ill A is reproduced, while saturability effects 
become important for larger amplitudes. Note in particu¬ 
lar that in the large-amplitude, strongly saturated, limit, 
f{x) —>■ 1 , so the system ( 1 ) becomes again effectively 
linear with just a frequency shift of the linear dispersion 
relation (3), A —)■ A -b 7 . 

We then obtain for the left-hand side in the general 
compacton condition (7): 

7 \A\^{i-i/m 

Jz(l + |Apy(l + |A/jp)’ 


which gives a second-degree equation for \A\^. Solving 
this and using Eq. ( 8 ) gives the following expression for 
the total compacton power, as function of the coupling 
ratio J and the rescaled nonlinearity/saturability param¬ 
eter Fs = 7 /J/: 


P+ = 


(i + 2/|J|2) 




± 






2-J2 


[jN -4|J|^ 


■ ( 16 ) 


It is important to stress that Fg and P act as two inde¬ 
pendent parameters for the saturable case, in contrast to 
the pure power-law case from Sec. Ill which only depends 
on the combination (q/Jjj^/'^P. Thus, a saturable non¬ 
linearity introduces an additional, qualitatively different, 
mechanism for compactification tuning via the saturabil¬ 
ity parameter Fg. 

We can note some general features of the expression 
(16): for given generic parameter values, there are ei¬ 
ther zero or two possible solutions with different power 
(but note of course that only solutions with P > 0 have 
physical meaning as compacton solutions). Also, for real 
J, the low power branch P_ and the high power branch 
P+ are positive within the same parameter space. Thus, 
an important difference to the cubic case is that the sat¬ 
urable nonlinearity promotes the appearance of two dif¬ 
ferent compaet modes for a given value of J, for the same 
existence regions. These two solutions possess the same 
profile (equal a), but different amplitude A. 

Bifurcation points, where two solutions coincide, ap¬ 
pear when the expression under the square root in (16) 
is zero, which happens when 


Ibif 


= rg 


(j^-2) (1±|J|) 

(It |J|) 


(17) 
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For a given r^, this gives a third-degree equation for the 
coupling ratios J where bifurcations appear. Note, in 
particular, that when Fg = — 2 , one bifurcation point ap¬ 
pears at J = 0 (with P —>• oo). The consequence is that, 
in addition to the two branches of solutions continuing 
the linear solution at J = two new branches appear 
for small J when Fg < —2 (see Fig. 6 (c), and lower 
left parts of (a) and (b)). Note also that at bifurcation 
points, the power is given by 

P± = |J| (i + 2/|J|2), 

and |y4p = | J|. 

From Fig. 6 we can see the limits of existence (non¬ 
black regions in top figures) of the compacton solutions 
vs. J and Fj. For F^ < —2, there are two windows of 
existence delimited by (white line) and Jun = 

For Fg > —2, there exists only one region of compactons 
delimited by and Jun = Two examples for the 
saturable compacton power (16) as a function of J for 
Fs = —2.5 (c) and Fg = 2.5 (d) are shown on the bottom 
of Fig. 6 , together with the corresponding expression for 
the cubic case (14), which only agrees with the saturable 
case in the small-power regimes, as it should. The re¬ 
gion in parameter space of J, where the lower saturable 
compacton branch P_ is well approximated by the cubic 
compactons, increases with increasing iFgl, as exempli¬ 
fied for Fg = ±10 in the insets of Figs. 7-8 (d). Note 
however that, in contrast to the cubic case, for each fixed 
Fg there is always an upper limit for | J| for compacton ex¬ 
istence along the branch with | J| > bifurcating from 
the linear compacton when Fg < 0. This limit, given by 
\Jmax\ = grows monotonically with [Fgl, and it 

follows from (17) that \Jmax\ ~ V—Fg for large negative 
Fg. Note also that, as Fg decreases towards large neg¬ 
ative values, the additional branch for small J born at 
Fg = —2 approaches the branch with |J| < 1 appear¬ 
ing for F < 0 in the cubic case (Sec. Ill A). However, in 
contrast to the cubic case also this branch has an up¬ 
per limit for | J| < 1, and it can be seen from (17) that 
\Jmax\ 2/|Fg| for large negative Fg. When Fg > 0, 

the condition to have P > 0 demands that compactons 
may only appear for | J| < v^, just as for the cubic case 
in Sec. Ill A. However, as seen in Fig. 6 , this solution 
branch now has a lower limit with | J| > I, and from (17) 
we now obtain \Jmin \ 1 ± 2 /rg for large positive Fg. 

Analogously to the studies for the cubic nonlinearity 
illustrated in Fig. 3, we have analyzed numerically the 
stability of the two saturable compacton solutions. The 
results are summarized in Fig. 6 (explicit examples for 
Fg = ±10 are shown in Figs. 7-8 (c)). When Fg < 0 (ef¬ 
fectively defocusing nonlinearity) and J > v^, the low 
power (P_) branch is, as in the defocusing cubic case, 
always stable, independent of the Fg-value. Decreasing 
Fg < 0, the high-power {P+) solutions destabilize for 
intermediate values of J, but generally regimes of sta¬ 
ble compactons remain for low and high J (the former 
corresponding to higher power and stronger saturation). 


0 0.01 





FIG. 6 : (Color online) Top (a) and (b): white and shaded 
(green) areas show the existence regimes for exact compactons 
with saturable nonlinearity where P > 0 from (16) (black re¬ 
gions imply no solutions). The non-vertical boundaries (white 
lines) show from (17). Additionally, we indicate the sta¬ 
bility g for (a) P- and (b) P+ solutions, where stable solutions 
are contained in the shaded (green) regions {g < 0 . 01 ), while 
white regions are unstable {g > 0.01). Bottom: two examples 
of power vs. coupling ratio J for Fg = —2.5 (c) and Fg = 2.5 
(d). Solid lines show the P- (lower red) and P+ (upper black) 
branches, while dashed blue line shows the corresponding ex¬ 
pression for the cubic case (14). 


Thus, for defocusing saturable nonlinearities two differ¬ 
ent stable compacton modes with different powers, having 
similar profiles (equal value of R) may exist for a given 
Fg. Fig. 6 also shows that both solutions appearing for 
J < 1 when Fg < —2 always are unstable. Similarly to 
the cubic case, this regime is not connected with flat- 
band linear modes. When Fg > 0 (effectively focusing 
nonlinearity), a similar scenario as in the focusing cubic 
case is observed for the compact mode continued from 
the linear limit at J = ^/2: it is stable for small powers 
but destabilizes through resonances with extended lat¬ 
tice modes for larger powers. However, similarly to the 























site n site n 

FIG. 7: (Color online) Diagrams for saturable nonlinear so¬ 
lutions with defocusing nonlinearity Fs = —10. (a) A versus 
J, (b) R versus J (c) g versus J, and (d) H versus J, with 
inset P versus J. Exact compactons are shown in Red (P-) 
and Black {P+) lines, while odd line localized solutions are 
shown in Blue and Orange for P = 1 and P — 10, respec¬ 
tively. Dashed lines correspond to odd tip solutions for same 
parameters. Profiles for different powers and J (marked by 
empty and filled squares and circles) are shown for (e)-(h) 
odd line and (i)-(l) odd tip modes. In (a) the linear (P —>■ 0, 
upper two bands) and saturated (P —> oo, lower two bands) 
linear band spectra are shown with shaded areas. 


defocusing saturable case, the compact mode generally 
restabilizes in the strongly saturated, high-power regime 
of the P+ branch, and thus we may also for the focus¬ 
ing case find regimes of two simultaneously stable com¬ 
pactons, although only in a narrow interval for J < y/2. 

To study the relation of these compact solutions with 
families of standard, non-compact, localized modes (cf. 
Fig. 4 for the cubic case), we first present in Fig. 7 
collected results for Fg = — 10 , i.e., an effectively de¬ 
focusing case with two different regions for analytical 
solutions. First, as discussed above, compact solutions 
are seen to bifurcate at J = from linear flat band 
modes with zero (P-) and infinite {P+) power, coincid¬ 
ing with the low and high amplitude (saturated) linear 
spectra (the latter being shifted to A -I- 7 due to sat¬ 


uration) as shown in Fig. 7(a). For J > ^/2, we see 
that both P_ and P+ compacton branches (but not the 
branches with J < 1) appear as compactifications of fam¬ 
ilies of standard localized modes with main localization 
at the central line site. Analogously to the cubic case in 
Sec. Ill A, we illustrate in Fig. 7 properties of numeri¬ 
cally obtained families of modes centered at line and tip 
sites, respectively, at two different fixed values of power: 
P = 1 and P = 10. Figs. 7(e)-(l) show profiles for low 
power solutions (empty squares and empty circles) and 
high power ones (black squares and black circles), for 
different values of J as well as for line and tip site cen¬ 
ters. Note that, while the low-power compactifying line- 
site centered mode (Figs. 7(e)-(f)) is always stable for 
the relevant regime in J, just as in the cubic case, the 
corresponding family of high-power modes for P = 10 
(Figs. 7(g)-(h)) has an instability window around the 
compactification point, and thus this compact mode is 
unstable. This instability window moves as the power 
is further increased, implying that, as mentioned above, 
also the high-power compacton branch stabilizes for a 
sufficiently high power (for Fg = — 10 , high-power com¬ 
pacton restabilization appears for P > 25.36, correspond¬ 
ing to J < 1.75). Note that also the non-compact modes 
with exponential tails are strongly localized (small R) for 
these parameter values, as they are residing far below the 
linear band. 

In Fig. 8, we illustrate in an analogous way the case 
with focusing nonlinearity, Fg = -1-10. As in the cubic 
case, neither the upper nor lower compacton branches 
then appear as compactifications of the family of fun¬ 
damental localized modes with constant sign and main 
localization at the central line site, which is seen to ap¬ 
pear simultaneously as a stable solution but with dif¬ 
ferent properties (Fig. 8 (g)). For small powers, when 
the compacton frequency is still in the linear gap, the 
scenario is similar to the cubic case illustrated in Fig. 5 
with a compactifying gap mode with a tail of oscillating 
signs (cf. Fig. 5 (f)). For slightly larger powers (P = 1 
in Fig. 8 ), when the compacton frequency has entered 
the upper linear band but still belongs to the P_ com¬ 
pacton branch, the compactifying family is not localized 
but has a non-decaying, oscillating tail (Fig. 8 (h)), which 
vanishes only at the compactification point. As seen in 
Fig. 8 (c), these solutions are generally unstable. Finally, 
when the power is large enough for the compacton fre¬ 
quency to reside above the linear spectrum (belonging 
then to the P+ branch for the case in Fig. 8 (a)), the 
compactifying family is again a localized solution with 
main localization at three sites, where only the central 
site has a different amplitude sign than the other sites 
(Fig. 8 (f)). From Fig. 8 (c), this solution is seen to be 
stable in parts of its existence regime (e.g., the mode in 
(f)), but for these parameter values, it is unstable at the 
point where it compactifies. However, for even larger P, 
at strong saturation, these solutions will also be stable 
at the compactification point, as illustrated in Fig. 6 (b). 
This restabilization of the eompactifying focusing mode 
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FIG. 8: (Color online) Diagrams for saturable nonlinear so¬ 
lutions with focusing nonlinearity, = 10. (a)-(d) show the 
same quantitites as Fig. 7 (a)-(d) with the same line coding 
for solutions. In addition, compactifying odd line localized 
solutions are shown in black dotdashed for P = 1 and green 
dotted for P = 10. Profiles at positions indicated by mark¬ 
ers: (e), (f) modes from the compactifying family at P = 10 
(localized modes above linear spectrum); (g) the fundamental 
(noncompactifying) stable odd line mode at P = 1; (h) mode 
from the compactifying family at P = 1 (non-localized mode 
inside upper linear band). Shaded areas in (a) represent dis¬ 
persion bands in the linear (P —>■ 0, lower two bands) and 
saturated (P —>■ oo, upper two bands) limits. 


above the linear spectrum at high power is another char¬ 
acteristic feature of the saturable system, which does not 
appear for the cubic case. 


V. FURTHER TUNING OPTIONS 


Let us now consider a generalization of Eq. (1), by 
adding a non-homogeneous distribution A„ of on-site en¬ 
ergies, and allowing also for anisotropies in the coupling 
coefficients Jn.m, 


iiln “b ^n'^n “b ^ ^ “b T/(|ttn| )^n — 0- (f^) 

m 


We were able to find three-site compact modes, having 
the same symmetric line-site centered structure as de¬ 
scribed in Sec. II, for a distribution of 


A„/J/ 


0 for n odd 

A for n even. 


(19) 


o— 


h 


■V 




-O 



FIG. 9: (Golor online) Sketch of the anisotropic geometry 
with pairwise alternating couplings and alternating on-site 
energies. 


This means, that the tips have an on-site energy of A, 
whereas the line sites have vanishing on-site energy (the 
linear dispersion relation and the corresponding linear 
compact mode in the isotropic sawtooth lattice with A = 
— 1 were illustrated in Ref. [8], Fig. 1 (f)). 

Considering also various ways to introduce coupling 
anisotropies in the sawtooth lattice, one possibility would 
be to have constant line coupling with alternating cou¬ 
plings to the tips. However, such a structure only allows 
for asymmetric linear compactons for some specific ratio 
of the two tip coupling constants, and no compact non¬ 
linear continuations were found for these modes for any 
distribution of on-site energies of the type (19) with iden¬ 
tical tip energies (nonlinear asymmetric compactons may 
exist under special conditions if the alternating coupling 
is also combined with different on-site energies for up- 
and down- pointing tip sites, but these conditions are 
more complex and will not be discussed further here). 
On the other hand, a nonlinear continuation of symmet¬ 
ric compact modes is possible, if instead a pairwise (real) 
alternating coupling Ji and J 2 is assumed in the tips, 
as sketched in Fig. 9. Within the line we conserve the 
coupling J; = 1. In this case, the four linear bands are 
given by 


A±,± = -^ + cos(2fc) ± 


( 20 ) 


]j{Jl + J 2 ) + (f" “ cos(2A:)^ ± 


with 


fpMk) = Jt + Jl - 2 JiV| + 4J2 J| cos2(2A) 


The geometry allows for two different types of com¬ 
pact modes, one that has the line site coupled to the tips 
with Ji, the other with J 2 . Without loss of generality we 
will consider only the latter, shown exemplary in Fig. 9. 
For a Kerr nonlinearity, the symmetric compact modes 
with amplitude ratio a. = —IjJi appear when the fol¬ 
lowing conditions for the effective nonlinearity (14) and 
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FIG. 10: (Color online) Compact mode Fa (solid), Aa 
( dashed), and stability parameter g (dotted) vs. J for A = 
—2.5 (-1.5,-0.5, 1.5) in (a)-(d). The shaded areas denote the 
band extensions. 


frequency are fulfilled. 


Ta{Ji,J2) 


Aa 


J2{Jt - 4) - JiA(Ji^ + 2) 

Ji - J! 

A J3 -f 2 J2 - Jfj2 

J! - Ji 


( 21 ) 


A. Alternating on-site energies 


We first turn to the isotropic case of Ji = J 2 = J, but 
A ^ 0. The linear dispersion relation (20) for this case 
is 


A±(fc) = Y-I- cos(2fc) ± 


A 


( 22 ) 


cos(2fc) —) -I-2 J2[l-I-cos(2fc)], 


whereas the Eqs. (21) simplify to 

(J2 -p 2) (J2 - A - 2) 


Ta = 


Aa = 


1-J2 

A J2 -h 2 - 
J^^l ■ 


(23) 


This gives a linear compact mode for JA,Zm = V^~+~A, 
thus there is only a continuation of the linear compact 
solution for A > —2. The bandgap vanishes as well for 
A = —2. However, nonlinear compact solutions exist 
also for A < —2. In Fig. 10 we show different scenarios. 
For A = —2.5 (see Fig. 10(a)), there is no flat band 
and the compact mode does not cross the bandgap. For 
A = —1.5 (see Fig. 10(b)), the flat band is located for 



FIG. 11: (Color online) Stability parameter g of the compact 
mode vs. Ji = J 2 = J and A (a), vs. Ji and J 2 with A = 0 
(b), only the black areas {g < 0.05) contain stable modes. 
N = 40. In (a) compactons are unstable also in the black 
stripe for small J when A > — 2, but have g < 0.05. The 
linear compacton condition is given by red lines, the ordinary 
sawtooth with homogeneous J 2 = Ji is marked with a gray 
line. 


J < 1, where the compacton frequency curve enters and 
crosses the bandgap with positive Ta. The sign of the 
discontinuity at J = 1 changes at A = — 1 and with 
it the location of the flat band, so for A = —0.5 and 
A = 1.5 (Fig. 10(c)-(d)) the crossing of the bandgap is 
found for J > 1. Increasing A further only opens the gap 
more, but there will be no further changes in regarding 
the bandgap. Therefore, the on-site energy difference A 
gives rise to the possibility of a direct engineering of the 
bandgap as well as the nonlinearity and frequency of the 
compacton. 

The stability for compact modes vs. J and A are shown 
with dotted lines for the specific examples in Figs. 10(a)- 
(d). The general stability scenario is explored in Fig. 
11 (a), where black areas denote stable or weakly unstable 
{g < 0.05) compactons and the degree of instability in¬ 
creases with brightness. As a guide to the eye, we denote 
the parameters of the linear compacton, A = J'^ — 2, 

with a red line. This line apparently separates two differ¬ 
ent regimes, with the main lobes of the instability located 
at opposite sides of J = 1. For a fixed A < —2 (Fig. 10 
(a)), there is always a stable regime for small J, corre¬ 
sponding to a compacton with main amplitude at the tip 
sites in a chain with focusing nonlinearity. For growing 
J, some weaker instabilities are passed for J < 1, and at 
J = 1 the large instability lobe appears as the compacton 
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now exists for defocusing nonlinearity. By comparing 
with the stability results for the corresponding solutions 
in the open trimer configuration (without a surrounding 
lattice) [26], it can be deduced that the instabilities on 
the focusing side (J < 1 ) are due to resonances with lat¬ 
tice oscillations (since this regime is always stable for the 
trimer [26]), while on the defocusing side (J > 1) it is 
an internal resonance in the compacton core which ap¬ 
pears also for the trimer. In the regime — 2 < A < — 1 
(Fig. 10 (b)), the main qualitative change appears in the 
small-J regime, where the nonlinearity now is defocus¬ 
ing for J < JA.iin s-nd the compacton is unstable. In 
fact, if follows from analogous considerations for the open 
trimer [26] that compactons with main localization on tip 
sites (small J) are always unstable through core instabil¬ 
ities for defocusing nonlinearities. For A > — 1 (Fig. 10 
(c)-(d)), the mode is always unstable in the defocusing 
regime 0 < J < 1. When A > —1, we also see an insta¬ 
bility in the other defocusing region for J > Ja./mj as a 
remains of the instability lobe at J > 1 for smaller A. 
For larger A, these instabilities disappear and the only 
significant instabilities observed for J > 1 are those in the 
focusing regime J < Jamu, which are analogous to those 
described for the homogeneous case A = 0 in Sec. Ill A. 
In general, for all A the compacton always approaches a 
stable one-site line-mode for large J, corresponding to a 
strong defocusing Kerr nonlinearity. 

One important conclusion from the above analysis 
is, that by introducing alternating on-site energies with 
A < —2 one may, for focusing nonlinearity, stabililize 
also compactons with main localization at the tip sites, 
which in the previous Sections were found to be always 
unstable for the homogeneous sawtooth lattice. 

B. Pairwise alternating couplings 

We finally turn briefly to the case of pairwise alternat¬ 
ing coupling as in Fig. 9, but considering for simplicity 
only vanishing on-site energies A = 0. As the linear sta¬ 
bility analysis vs. Ji and J 2 for this family shows (see Fig. 
11 (b)), for small external coupling to the array Ji < 1 
(corresponding to defocusing nonlinearity, F < 0 accord¬ 
ing to ( 21 )) nearly all modes (having their main ampli¬ 
tudes on tip sites) are unstable, only a very restricted 
zone for small J 2 < Ji is stable. On the other hand, 
if Ji > 1, there are as before many regimes of stable 
modes. In the focusing regime, 1 < Ji < y/2 (see (21)), 
the stability scenario is qualitatively similar as for the 
case J 2 = Ji described in Sec. Ill A; changing the internal 
coupling J 2 mainly changes the internal compacton oscil¬ 
lation frequencies and thereby also the exact location of 
the instability threshold resulting from their resonances 
with the surrounding chain, but a stable regime generally 
remains for small nonlinearity. In the defocusing regime 
Ji > the compacton having main amplitude on the 


line-site always remains stable for J 2 > Ji, but it may 
destabilize in certain regimes for J 2 < Ji, i.e., when the 
internal coupling is weaker than the external. Note in 
particular the instability window around J 2 = v^, cor¬ 
responding to a resonance with a linear compacton with 
internal coupling Ji. 


VI. CONCLUSIONS 

In this work we showed the existence of compact modes 
for generic realizations of the sawtooth lattice, consider¬ 
ing power law and saturable nonlinearities. These modes 
are compact nonlinear continuations of the linear flat- 
band modes found in the sawtooth model. For a wide 
window in parameter space we showed the stability of 
these modes for both focusing and defocusing nonlinear¬ 
ities. Furthermore, we explored the influence of inho¬ 
mogeneous on-site energies and pairwise alternating cou¬ 
pling on the band-structure and stability of the compact 
modes. We observed, that negative on-site energies or 
high anisotropic pairwise coupling increase the stability 
regimes. 

Moreover, we have confirmed the existence of com¬ 
pacton solutions for the saturable nonlinearity, finding 
stable regimes for two types of analytical solutions. As 
they are very compact solutions, they could be easily 
excited in an experiment, for low level of power and 
for different sawtooth geometries governed by the cou¬ 
pling ratio J. The possibility to excite very localized 
nonlinear solutions using low level of power is an im¬ 
portant property of this particular lattice, in contrast 
with nonlinear conventional lattices [14]. We also found 
regimes with possibility to excite two, simultaneously sta¬ 
ble, compactons with similar profiles, for small and large 
level of power. This is a very interesting property of 
the sawtooth lattice; in conventional systems, stable so¬ 
lutions for different level of power exist but with different 
profile structure (e.g., different participation number). 
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